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Measuring the prompt atmospheric neutrino flux with down-going muons in neutrino 

telescopes 
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In the TeV energy region and above, the uncertainty in the level of prompt atmospheric neutrinos 
would limit the search for diffuse astrophysical neutrinos. We suggest that neutrino telescopes may 
provide an empirical determination of the flux of prompt atmospheric electron and muon neutrinos 
by measuring the flux of prompt down-going muons. Our suggestion is based on the consideration 
that prompt neutrino and prompt muon fluxes at sea level are almost identical. 



Atmospheric neutrinos and muons, i.e. neutrinos and 
muons produced in the atmosphere by cosmic ray inter- 
actions, are the most important source of background 
for present and future high-energy neutrino telescopes, 
which are expected to open a new window in astronomy 
by detecting neutrinos from astrophysical sources [1] . (In 
this Letter, 'muons' includes 'antimuons' and 'neutrinos' 
includes 'antineutrinos'.) 

In their current design, neutrino telescopes consist of 
large arrays of phototubes located under water or ice. 
They detect high-energy neutrinos through the charged 
particles these produce in the water or ice inside or 
around the instrumented array. 

Atmospheric muons can reach the detector only from 
above, because the range of muons in Earth is only a few 
kilometers. Atmospheric muons are therefore only down- 
going. Their flux is typically so high that the region of 
sky accessible to even very deep neutrino telescopes is 
only the hemisphere below the horizon. 

Atmospheric neutrinos can instead reach the detector 
from all directions. Hence they are an irreducible back- 
ground for diffuse astrophysical neutrino fluxes. It is 
therefore very important to evaluate their intensity with 
reasonable accuracy. 

At GeV energies the atmospheric muon and neutrino 
fluxes are dominated by 'conventional' sources, i.e. de- 
cays of relatively long-lived particles such tt and K 
mesons. With increasing energy, the probability in- 
creases that such particles interact in the atmosphere be- 
fore decaying. This implies that even a small fraction of 
short-lived particles can give the dominant contribution 
to high energy muon and neutrino fluxes. These 'prompt' 
muons and neutrinos arise through semi-leptonic decays 
of hadrons containing heavy quarks, most notably charm. 

Estimates of the magnitude of the prompt atmospheric 
fluxes differ by almost 2 orders of magnitude. Fig. 1 
shows a compilation of prompt muon fluxes at sea level. 
Prompt neutrino fluxes are essentially identical, while 
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FIG. 1. Vertical atmospheric muon (and neutrino) fluxes. 
Muon and neutrino conventional fluxes from [2] (below 10^ 
GeV) and [3] (above) (line marked conven. and dashed lines). 
Muon prompt fluxes from: two charm production models in 
[4] (ZHVa, ZHVe); empirical model in [5] (RVS); quark-gluon 
string model and recombination quark-parton model in [6] 
(QGSM, RQPM); perturbative QCD in [7] (PRS band), [8,9] 
(GGV band), and [3] (TIG). Also indicated are the experi- 
mental bounds on prompt muons from LVD [10] and AKENO 
[11]. 



conventional neutrino fluxes are lower by one (i/^) or two 
(I'e) orders of magnitude. The crossing from conventional 
to prompt muon fluxes happens between 40 TeV and 
3 PeV, while the analogous crossing for muon-neutrinos 
happens at lower energies, between 20 TeV and 800 TeV. 
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The uncertainty in the intensity of conventional atmo- 
spheric neutrinos and muons is thought to be approxi- 
mately 30% at present, but could decrease to about 10% 
with coming improvements [12]. At about 1 TcV, the 
contribution of prompt neutrinos taking into account the 
LVD bound could be as high as 10% of the conventional 
neutrino flux. 

Between 1 TeV and 100 TeV, prompt neutrinos be- 
come the biggest source of uncertainty in the atmospheric 
neutrino flux. So the level of prompt neutrinos is a po- 
tential problem which would limit the search for diffuse 
astrophysical neutrinos at energies of about 1 TeV, much 
smaller than the energies where they become dominant 
(see for example [12]). 

Here we suggest a way to overcome the theoretical un- 
certainty in the magnitude of the prompt electron and 
muon neutrino fluxes by deriving their intensity from a 
measurement of the down-going prompt muon flux. Our 
suggestion is based on the observation that, due to the 
charmed particle decay kinematics and the same branch- 
ing ratios for the semi-leptonic decays into efg and fiv^, 
the prompt electron and muon neutrino fluxes and the 
prompt muon flux are essentially the same at sea level 
[3,7-9,13]. This result is independent of the charm pro- 
duction model. 

We want to stress that we are suggesting to use down- 
going prompt muons and not up-going neutrino-induced 
muons whose flux is orders of magnitude smaller. While 
an important contribution to up-going muons is expected 
from astrophysical neutrinos, no astrophysical signal is 
expected in down-going atmospheric muons. 

Moreover, prompt muons are much easier to detect 
than prompt neutrinos, since the latter have to convert 
to a charged particle within the effective volume of the 
detector. In fact, the flux of upcoming muons induced by 
muon neutrinos at 1 TeV is about 10"'' of the neutrino 
flux at sea level (using charge-current cross sections in 
[14] and muon ranges in [15]). On the other hand, the 
flux of down-going muons at a slant depth of about 3 km 
w.c. is only a fraction 0.4 0.6 of the muon flux at sea 
level (the exact suppression factor depending on depth, 
energy spectra and zenith angle [15]). Thus close to 1 
TeV, taking into account that the conventional neutrino 
fluxes at sea level are about 10% of the muon fluxes (see 
Fig. 4 below), for each up-going neutrino-induced muon 
there are 4 6 lO'' down-going muons. Of these, as much 
as a few percent may be prompt. 

From Fig. 1, the vertical conventional muon flux above 
1 TeV is approximately 10^" muons/km^/yr/sr at sea 
level. This implies roughly 10^ down-going events per 
year in a kilometer-size detector at a depth of about 3 
km. Thus, to extract a 1% fraction of prompt muons at 
1 TeV, it would sufSce to record 1 out of 10^ down-going 
events per bin in the sky for a year (fewer events would 
need to be recorded at higher energies). 

For what we suggest, it is important to separate the 



prompt muons from the conventional ones for two rea- 
sons: (1) the conventional neutrino fluxes are small frac- 
tions, less than 10%, of the conventional muon flux, and 
our method of using the ratio of neutrino to muon fluxes 
would become less straightforward; and (2), as a con- 
sequence of the previous reason, the ratio of neutrino to 
muon fluxes depends on the crossing energy between con- 
ventional and prompt fluxes, and so on the large uncer- 
tainty on the absolute value of the prompt fluxes, making 
our method inapplicable. 

There are ways of separating the prompt muons from 
the conventional ones in underwater or under-ice detec- 
tors, such as the different zenith angle dependence of the 
prompt and conventional fluxes; the different depth de- 
pendence at a given zenith angle; and the different spec- 
tral shape at a given depth and zenith angle (see e.g. 
Ref. [16,17]) 

In a series of papers [8,9,13] (called GGVl, GGV2 

and GGV3 from now on), we studied the prompt lep- 
ton fluxes using a model for charm production in the at- 
mosphere based on Quantum Chromo-Dynamics (QCD), 
the theoretically preferred model. We used a next- 
to-leading order perturbative QCD (NLO pQCD) cal- 
culation of charm production, as implemented in the 
Mangano-Nason-Ridolfl program [18] calibrated at low 
energies, followed by a full simulation of particle cascades 
in the atmosphere generated with PYTHIA routines [19]. 

In our first paper (GGVl), we tried different modes of 
cascade generation, different options allowed by PYTHIA 
in the various stages of parton showering, hadronization, 
interactions and decays, etc., noticing changes of at most 
10% in the final results. We decided to use what we 
called our 'single' mode simulation, with showering, in- 
dependent fragmentation, interactions and semileptonic 
decays according to Ref. [3]. In our 'single' mode we 
enter only one c quark in the particle list of PYTHIA, 
and we miiltiply the result by a factor of 2 to account 
for the initial c quark. PYTHIA performs the shower- 
ing, standard independent fragmentation, and follows all 
the interactions and decays using default parameters and 
options. 

In GGVl we found that the NLO pQCD approach pro- 
dTices fluxes in the bulk of older predictions (not based 
on pQCD) as well as of a pQCD semianalytical analysis 
[7] . We also explained the reason of the low fluxes of the 
model of Ref. [3] , which were due to the chosen extrapo- 
lation of the gluon partonic distribution function (PDF) 
at small momentum fractions x. 

In GGV2, we considered four sets of PDF's: MRS Rl- 
R2 [20], CTEQ 4M [21] and MRST [22]. Besides the 
choice of the PDF set, our procedure has the freedom 
to choose reasonable values of the charm mass m,c, the 
factorization scale /j-f, and the renormalization scale /z_r, 
so as to fit the experimental data. In GGVl and GGV2 
we made the standard choice [18,23] of fxp = 2mT, Mil = 
rriT, where rriT — \Jpt + is the transverse mass. The 
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FIG. 2. Ratio of prompt neutrino to prompt muon fluxes 
as a function of lepton energy using the MRST PDF 
with A = (sohd line), 0.1 (dotted), 0.2 (short-dashed), 
0.3 (long-dashed), 0.4 (short-dashed dotted), and 0.5 
(long-dashed dotted). 



values of the charm mass were taken slightly clifFcrcnt 
for each PDF sot, namely: = 1.185 GeV for MRS 
Rl. m,c = 1.310 GeV for MRS R2, = 1.270 GeV for 
CTEQ 4M, and = 1.250 GeV for MRST. 

Due to the steep decrease with increasing energy of 
the incoming flux of cosmic rays, only the most ener- 
getic charm quarks produced count and those come from 
the interaction of projectile partons carrying a large frac- 
tion of the incoming nucleon momentum. Thus, the 
characteristic x of the projectile parton, xi, is large, 
Xi ~ 0(10^^). We can then immediately understand 
that very small parton momentum fractions are involved 
in pQCD charm production as follows. Typical partonic 
center of mass energies Vs are close to the cc threshold 
2mc — 2 GeV, (since the differential cc production cross 
section decreases with increasing s) while the total center 
of mass energy squared is s = 2mNE (with ttiat — 1 GeV 
the nucleon mass, and E the energy per nucleon of the in- 
coming cosmic ray). Calling X2 the momentum fraction 
of the target parton in the nucleus of the atmosphere, 
we have xiX2 = s/s = 4m^/(2mAri?) ~ GeV /E. Hence 
Xi ~ O{GeY/0.lE) ~ 0{GeV/Ei), where Ei ~ O.IE is 
the dominant muon or neutrino energy. 

In GGV2, we analyzed in detail the dependence of the 
fluxes on the extrapolation of the gluon PDF at small x, 
which, according to theoretical models, is assumed to be 
a power law with exponent A, xg{x) ~ x~^, with A in the 
range 0-0.5. Particle physics experiments are yet unable 
to determine the value of A at a; < 10~^. We found that 
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FIG. 3. As Fig. 2, but for CTEQ 4M and A = and 0.5. 



the choice of different values of A at a; < 10 ^ leads to a 
wide range of final prompt fluxes at energies above 10^ 
GeV. 

Due to this result, in GGV2 and GGV3 we suggested 
the possibility of measuring A through the atmospheric 
muon fluxes at energies above 10^ GeV, using not the ab- 
solute fluxes, because of their large theoretical error, but 
rather their spectral index (i.e. the "slope" of the flux). 
In particular, in GGV3 we proposed to use the slope of 
the flux of down- going prompt muons, and presented an 
overall error analysis of the model we used. 

Here again we are suggesting to use down-going 
prompt muons, at energies iJ^ > 1 TeV where prompt 
muons can be separated from conventional ones [17], this 
time to measure the flux of prompt electron and muon 
atmospheric neutrinos at sea level. We find that the ratio 
of prompt neutrino to prompt muon fluxes is about 1.1, 
constant with energy to within 10%, and almost indepen- 
dent of the choice of PDF and charm production param- 
eters. This is shown in Figs. 2 and 3 for the MRST and 
CTEQ 4M PDF's with a range of A values from to 0.5. 
We do not show results obtained with the MRS Rl and 
R2 PDF's because they are similar. We expect that other 
models of charm production in the atmosphere, even not 
based on perturbative QCD, will lead to a similar ratio, 
because this ratio depends essentially only on the decay 
properties of the charmed hadrons. 

To complete the discussion, in Fig. 4 we plot the 
neutrino-over- muon ratio of the sum of conventional plus 
prompt lepton fluxes as a function of lepton energy. In 
this figure, we use the conventional fluxes in Ref. [3] and 
the prompt fluxes of GGV2 with the MRST PDF and 
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FIG. 4. Total neutrino-over-muon ratio as a function of 
lepton energy. Vertical marks denote the crossing energy from 
conventional to prompt muons. 



A = (thick lines) and 0.5 (thin hnes). For each A, 
the crossing energy from conventional to prompt muons, 
from Fig. 2 in GGV2, is marked with a vertical stroke. 

Wc have suggested a way to overcome a potential prob- 
lem which would limit the search for diffuse astrophysical 
neutrinos in underwater or under-ice neutrino telescopes, 
namely the theoretical uncertainty of about two orders 
of magnitude in the intensity of the prompt atmospheric 
neutrino fluxes. Concretely, we have suggested to deter- 
mine their intensity from a measurement of the down- 
going prompt muon flux at sea level, whose intensity is 
the same to within 10% or better. 
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